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Aims. We probe the evolution of the fine-structure constant a with cosmic time. 



> ABSTRACT 

(N 

, Methods. Accurate positions of the Fell lines A1608, A2382, and A2600 are measured in the z a b s = 1.84 absorption system 
^—^ from a high-resolution (FWHM ~ 3.8 km s~ ) and high signal-to-noise (S/N > 100) spectrum of the quasar Q 1101-264 
(zcm = 2.15, V = 16.0), integrated for 15.4 hours. The Single Ion Differential a Measurement (SIDAM) procedure and the Ax 2 
, method are used to set constraints on Aa/a . 
Results. We have found a relative radial velocity shift between the A1608 and AA2382, 2600 lines of Ad = —180 ± 85 m s" 1 
(both random and systematic errors are included), which, if real, would correspond to Aa/a = (5.4 ± 2.5) x 10~ 6 (la C.L.). 
Considering the strong implications of a such variability, additional observations with comparable accuracy at redshift z ~ 1.8 
^ ' are required to confirm this result. 

' Key words. Cosmology: observations - Line: profiles - Quasars: absorption lines - Quasars: individual: Q 1101-264 
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1. Introduction scale distribution of galaxies, and the integrated Sachs- 



m , , . , , , , . , r . „ , , ^ , Wolfe effect (for a review, see Copeland et al. 2006). 

The current cosmological model (ACDM, for A Cold Dark v ' 

Matter) assumes that our universe is dominated by dark The observational data on SNe la suggest that the 
matter and dark energy. This model is based on a series of kinematics of the expansion of the universe, decelerating 
independent observations such as the luminosity distances at red shift z > 1, has changed to the accelerating phase at 
of high-redshift type la supernovae (SNe la), the cosmic lower z - This late-time acceleration of the universe is re- 
microwave background temperature anisotropy, the large S arded as evidence for a non-zero cosmological constant A 
or for the existence of dark energy with negative pressure. 
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* Ine nature ol tne dark energy is still unknown. At 
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* „ , , .. , , . ,rTm T ^ ^1 the moment, a variety ol models is being considered, from 

Based on observations performed at the VL1 Kueyen tele- . J; 

/-i^c/-. t-> i <^u-i \ j.v. t^c/^ m nTC A scalar fields with dificrcnt potentials to extra dimensions 

scope (ESO, Paranal, Chile), the ESO programme No. 076. A- K 

Q4g 3 (D > 4) of physical space-time (e.g., Dvali et al. 2000). In 

** On leave from the Ioffe Physico-Technical Institute, St. man y models, the cosmological evolution of dark energy 

Petersburg, Russia is accompanied by variations in coupling constants, the 

*** On leave from the Ioffe Physico-Technical Institute, St. fine-structure constant a and the electron-to-proton mass 

Petersburg, Russia ratio. This means that the evolution of dark energy — if it 
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indeed occurs — can be traced through the measurements 
of the coupling constants changes in course of cosmic time 
(e.g., Avelino et al. 2006). Theories predict different be- 
havior of dark energy, from slow- rolling to oscillating (e.g., 
Liddle & Lyth 2000; Fujii & Maeda 2003). It is clear that 
in order to study the dynamics of dark energy the fun- 
damental constants must be measured at every redshift z 
with highest accuracy. 

The value of Aa/a= (a z — a) /a at a redshift z can 
be estimated from the radial velocity shifts of the line po- 
sitions of different ions. Theoretical basics of the method 
called 'Many Multiplet (MM) method' are given in Webb 
et al. (1999) and Dzuba et al. (2002). Recently obtained 
results are as follows. 

Murphy et al. (2004) claim that Aa/a= -5.7 ± 1.1 
ppm (ppm stands for parts per million, 10 -6 ) in the red- 
shift range 0.2 < z < 4.2, i.e. a was lower in the past. This 
result was produced by averaging over 143 absorption sys- 
tems detected in HIRES/Keck telescope spectra of QSOs. 
However, Chand et al. (2004) report Aa/a= -0.6 ± 0.6 
ppm for 0.4 < z < 2.3 from the analysis of 23 absorption 
systems identified in the UVES/VLT spectra of QSOs. It 
should be emphasized that 

- both the Murphy et al. and Chand et al. estimates are 
obtained by averaging over many absorption systems 
with different rcdshifts, i.e. the accuracy of Aa/a at 
individual redshifts is considerably lower than the re- 
ported values; 

Mg II and Fe II ions — the main elements of the most 
accurate low-z sample by Murphy et al. — cannot be 
formed in the same volume of the intervening metal ab- 
sorbers as shown from spectral observations of GRBs 
(Hao et al. 2006; Porciani et al. 2007). 

Recent laboratory restriction on the rate of fractional 
variation of the fine-structure constant of (—2.7 ± 2.6) x 
10~ 15 yr _1 (Cingoz et al. 2007), which is based on highly 
sensitive to temporal variation of a transition frequencies 
in two isotopes of dysprosium (Dy), indicates no signif- 
icant variation at present time (z = 0). Revision of the 
fission products of a natural reactor in Oklo (z = 0.4) 
also indicates no variation at the level of 1.2 x 10 -17 yr -1 
(Gould et al. 2006). 

In order to improve the accuracy of individual 
Aa/a values measured from quasar absorption-line spec- 
tra, we developed a modification of the MM method 
called 'Single Ion Differential a Measurement, SIDAM' 
(Levshakov 2004; Quast et al. 2004; Levshakov et al. 
2005, hereafter Paper I; Levshakov et al. 2006, hereafter 
Paper II). Based on the analysis of only one element 
(singly charged iron, Fen), it is free from the influence 
of small Doppler shifts between different elements caused 
by ionization inhomogcneities within the absorber, and by 
unknown isotopic abundances. It is also free from different 



line profile distortions due to finite correlation length ef- 
fects in turbulent velocity fields which are sensitive to the 
ratio of the rms turbulent velocity to the thermal width of 
the absorbing ions (Levshakov & Kegel 1997; Levshakov 
et al. 1997). 

Using SIDAM, we have reached the accuracy of 
Aa/a for one absorption system at the level of <~ 10~ 6 , 

1. e. as high as the accuracy of the mentioned above val- 
ues averaged over many absorption systems. Our previous 
results give the following constraints on Aa / a . 

At z abs = 1.15 towards HE 0515-4414, Aa/a = -0.4± 
1.9 ppm at z a bs = 1-15 (Quast et al. 2004). Later on, we 
improved the accuracy of this estimation and obtained 
Aa/a= -0.07 ± 0.84 ppm (Paper II). This utmost ac- 
curacy was achieved due to unique favorable conditions 
such as the brightness of the quasar and the strength of 
all Fen lines including A1608 A. In both cases the same 
archive data were used: FWHM ~ 5.6 km s _1 in the blue 
and ~ 4.8 — 5.5 km s _1 in the red. A consistent restriction 
on Aa/a in this system was obtained independently by 
Chand et al. (2006), Aa/a = 0.5 ± 2.4 ppm, who used the 
high resolution pressure and temperature stabilized spec- 
trograph HARPS mounted on the ESO 3.6 m telescope at 
the La Silla observatory. 

Further, using the VLT/UVES archive spectra of Q 
1101-264 with the resolution FWHM ~ 6.0 km s" 1 in the 
blue (A - 4500 A), and ~ 5.4 km s" 1 in the red (A - 7000 
A) we obtained Aa/a — 2.4 ± 3.8 ppm at z a b s = 1.839 
(Paper I). 

Recently this quasar was re-observed with higher spec- 
tral resolution. In the present paper we describe the new 
measurement of Aa/a at z a b s = 1.839 which allowed us to 
enhance the accuracy of Aa/a by a factor 1.5. 

2. Observations and Data Reduction 

The observations of Q 1101-264 were acquired with the 
VLT UV- Visual Echelle Spectrograph (UVES) on the 
nights of 2006 February 21, 22, and 23. The selected expo- 
sures are listed in Table 1 . The spectra were recorded with 
a dichroic filter which allows to work with the blue and 
red UVES arms simultaneously as with two independent 
spectrographs, and the CCDs were read out unbinned. 
The standard settings at central wavelengths A437 nm and 
A860 nm were used for the blue and red arms respectively 
(marked by symbols 'b' and 'r' in Col. 1 of Table 1). The 
slit was oriented along the parallactic angle and the width 
set at 0.5 arcsec, with DIMM values ranging from 0.5 to 
1.1 arcsec during the run. 

Calibration exposures were taken immediately after 
scientific exposures to minimize the influence of chang- 
ing ambient weather conditions which may cause different 
velocity offsets in the lamp and QSO spectra if they were 
not obtained closely in time (Paper II). Variations in tern- 
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Table 1. UVES/VLT observations of Q 1101-264 (ESO 
programme 076.A-0463) 



Exp. 


Date 


Time 


Exp. 


S/N per 


No. 


y-m-d 


h:m:s 


s 


pixel 


(1) 


(2) 


(3) 


(4) 


(5) 


lb 


2006-02-21 


07:39:26 


6400 


30 


lr 


2006-02-21 


07:39:22 


5400 


32,34 


2b 


2006-02-22 


02:28:57 


5400 


31 


2r 


2006-02-22 


02:28:52 


5400 


39,37 


3b 


2006-02-22 


04:04:47 


4953 


30 


3r 


2006-02-22 


04:04:43 


4959 


44,44 


4b 


2006-02-22 


05:55:21 


5400 


34 


4r 


2006-02-22 


05:55:17 


5400 


47,49 


5b 


2006-02-22 


07:30:28 


6708 


38 


5r 


2006-02-22 


07:30:25 


6712 


52,48 


6b 


2006-02-23 


01:29:26 


6300 


26 


6r 


2006-02-23 


01:29:25 


6300 


36,38 


7b 


2006-02-23 


03:19:23 


7200 


35 


7r 


2006-02-23 


03:19:22 


7200 


46,48 


8b 


2006-02-23 


05:45:15 


5400 


37 


8r 


2006-02-23 


05:45:10 


5400 


41,46 


9b 


2006-02-23 


07:19:40 


7500 


40 


9r 


2006-02-23 


07:19:34 


7500 


52,53 



Note: in Col. 1, 'b' and V denote the blue-arm (echelle 
order 102) and red-arm (echelle orders 90 and 83) 
exposures; calibration lamp exposures were taken 
immediately after each science exposure; the pixel 
sizes close to the Fell lines, where the S/N has been 
computed, are 25 mA (A1608), 30 mA (A2382), and 
27 mA (A2600). 

perature can induce color effects in the spectra since they 
act differently on the different cross disperscrs. 

To control possible systematics in radial velocities, 
the thermal and pressure stabilities were monitored. 
Variations of the ambient temperature during science ex- 
posures (the difference AT between the beginning and 
the end) were always less than 0.1 K. This is also the 
maximum variation occurring between the science expo- 
sure and the calibration lamp. The pressure variation AP 
reached 0.6 mbar in one case (No. 4, Table 1), in all other 
cases the variations of AP did not exceed 0.3 mbar. The 
estimations of Kaufer et al. (2004) arc of 50 m s _1 for 
AT = 0.3 K and AP = 1 mbar, so the observed differ- 
ences should not affect significantly the radial velocities. 

From the blue frame we used only order 102, and from 
the red — orders 90 and 83 (Col. 5), where Fell lines with, 
respectively, A = 1608, 2382 and 2600 A are observed. 
The selected lines are located close to the central regions 
of the corresponding echelle orders. This minimizes possi- 
ble distortions of the line profiles caused by the decreas- 
ing spectral sensitivity at the edges of echelle orders. The 
signal-to-noise ratios (S/N) per pixel (25 mA, 30 mA and 
27 mA at the positions of A1608, A2382 and A2600, respec- 



tively) for the individual exposures and individual orders 
are listed in Col. 6. Cols. 2, 3 and 4 give, respectively, the 
date of observation, the beginning of the exposure, and 
the total exposure time. 



The spectra were reduced using the UVES context 
within the ESO MIDAS data reduction software. Master 
fiat-field and master bias frames were obtained by com- 
puting the median of five flat-field and bias frames taken 
in the afternoon. The master bias was subtracted from 
each science frame; subsequently the background was sub- 
tracted by fitting a two-dimensional spline to points in the 
inter-order region, where the background was estimated as 
the median number of counts over a box 51 pixels wide in 
the dispersion direction and 5 pixels wide in the cross- 
dispersion direction. The background subtracted science 
frame was divided pixel-by-pixel by the normalized mas- 
ter flat-field. The spectrum was extracted using an aver- 
age extraction and an optimal extraction. The S/N ratio 
proved to be higher in the average extraction, therefore 
this was used for the analysis. The wavelength calibration 
was derived from the spectrum of a ThAr hollow cath- 
ode lamp. We used the Los Alamos table of ThAr lines 
(Palmer & Engleman 1983; de Cuyper & Hensberge 1988) 
and the ID option, by which a 4th order polynomial is fit 
to the measured positions of arc lines in each echelle order. 
The wavelength calibration was applied to the extracted 
science spectrum both on a pixel by pixel basis, and re- 
binning the extracted spectrum to a constant wavelength 
step. We did not observe any systematic shift in the line 
centroids introduced by the rebinning, therefore we used 
the rebinned spectra, which can be readily co-added. We 
kept a very small tolerance in the fitting in order to reject 
several faint lines with lower accuracy (about 150 m s _1 ) 
still preserving somewhat 20 lines for each order, normally 
strong lines with about 15 m s _1 accuracy. Finally, as a 
cross-check, we performed an independent reduction of the 
spectra using a 2D wavelength calibration, by which a 2D 
polynomial dispersion relation of the form /(m, A), where 
m is the absolute order number and A is the wavelength, is 
fitted to the entire echellogram of the ThAr lamp (Lopez 
et al. 2005). The residuals of the calibrations in both ap- 
proaches were rather small corresponding to the rms errors 
of the wavelength scale calibration of about 20 m s _1 . The 
observed wavelength scale of each spectrum was trans- 
formed into vacuum, heliocentric wavelength scale (Edlcn 
1966). The resulting spectral resolution as measured from 
the ThAr emission lines was FWHM ~ 3.95 km s _1 in 
the blue (A ~ 4500 A) and ~ 3.75 km s _1 in the red 
(A ~ 7000 A). Thus, the resolution of these new spectra 
is significantly higher than in the previous observations (6 
and 5.4 km s _1 respectively, see Paper I). 
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Fig. 1. Combined absorption-line spectra of Fen asso- 
ciated with the z = 1.84 damped Lya system towards 
Q 1101-264 (normalized intensities are shown by his- 
tograms). The zero radial velocity is fixed at z a ^ s = 
1.838911. The synthetic profiles are over-plotted by the 
smooth curves. The normalized residuals, {Tf al -J=° hs 
are shown by dots (horizontal dotted lines restrict the 
la errors). The dotted vertical lines mark positions of 
the sub-components. Bold horizontal lines mark pixels 
included in the optimization procedure. The ranges at 
v < -50 km s^ 1 and at v ~ -30 km s^ 1 in the Fen A2600 
profile are blended with weak telluric lines. The normal- 
ized xl = 0-901 (i/ = 257). 

3. Estimate of Aa/a 

The spectroscopic measurability of Aa/a is based on the 
fact that the energy of each line transition depends in- 
dividually on a change in a (e.g., Webb et al. 1999). It 
means that the relative change of the frequency luq due 
to varying a is proportional to the so-called sensitivity 
coefficient Q — q/u)Q (Paper I). The q- values for the reso- 
nance UV transitions in Fe II are calculated by Dzuba et 
al. (2002). The corresponding sensitivity coefficients are 
listed in Table 5 in Paper II. The rest frame wavelengths 
for some Fen transitions have been re- measured recently 
with higher accuracy by Aldenius et al. (2006). For our 
present work we need two of them: 2382.7641 ± 0.0001 A, 
and 2600.1722 ±0.0001 A. The line from the 8th UV mul- 
tiplet 1608.45080 ± 0.00008 A in Table 5 is registered to 



Fig. 2. Same as Fig. 1 but for the 11-component model 
fitted to the central part of the Fe n lines. The normalized 
xl = 0.720(^=188). 

the calibration scale defined in Norlen (1973), whereas the 
lines A2382 and A2600 from Aldenius et al. are calibrated 
in the scale of Whaling et al. (1995). The two scales differ 
by a factor of about 7 x 10~ 8 in frequency space. Thus, 
a corrected wavelength of the A1608 line is 1608.45069 A. 
Other atomic data are taken as they are in Table 5. 

The value of Aa/a itself depends on a proper inter- 
pretation of measured relative radial velocity shifts, Av, 
between lines with different sensitivity coefficients, AQ. It 
can be shown that in linear approximation (|Aa/a| <1) 
these quantities are related as (Paper II): 

Aa __ (v 2 - vi) _ Av 

a ~ 2c(Q 1 -Q 2 ) ~ 2cAQ ' 1 ' 

Here index '1' is assigned to the line A1608, while index 
'2' marks one of the other Fen lines (A2382 or A2600). 

It is important that AQ must be large enough to pro- 
vide smaller uncertainty in Aa/a. With Qi = —0.021 
and Q 2 = 0.035 for both A2382 and A2600 lines, we 
have AQ = -0.056. Note that |AQ| = 0.056 is almost 
two times larger than |AQ| ~ 0.03 which is used in the 
standard MM method when Mgn AA2796,2803 and Fen 
AA2344, 2374, 2382, 2586, 2600 transitions are compared^. 

1 For Mgn, Q 2 796 = 0.00 59, and Q 2 803 = 0.0034. For Fell, 
Q2344 = 0.028, Q2374 = 0.038, and Q 25 86 = 0.039. 
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Fig. 3. Same as Fig. 2 but for the 10-component model. 
The normalized xl = 0.871 (y = 191). 

It is also important that the Fell transitions at A1608 and 
AA2382, 2600 used in the present study have sensitivity co- 
efficients of different signs which helps to distinguish the 
influence of hidden blends on the line position measure- 
ments (Levshakov & D'Odorico 1995). 

The individual vacuum-barycentric spectra were re- 
sampled to equidistant wavelength grids using smallest 
pixel sizes in the blue and red exposures and applying lin- 
ear interpolation. The resulting normalized and co-added 
spectra from 9 exposures (Table 1) are shown in Fig. 1. A 
weighted mean signal-to-noise ratio per pixel of S/N — 100 
and 120 was achieved in the final spectrum at A ~ 4566 A 
and A ~ 6765, 7382 A, respectively. Two independent data 
reduction procedures resulted in almost identical co-added 
spectra. In contrast to the archive data used in Paper I 
which we failed to fit self-consistently, this time we found 
a good fit of Fen profiles to a model, the normalized x 2 
per degree of freedom equals xl — 0.901 (y = 257). The 
computational procedure was the same as in Paper II. The 
optimal solution was found for the 16-component model 
shown in Fig. 1 by the smooth curves (the positions of the 
sub-components are marked by dotted lines). The model 
parameters are given in Table 2. Since the Q values for 
A2382 and A2600 are equal, their relative velocity shift 
characterizes the goodness of wavelength calibrations of 
the echelle orders 90 and 83. The revealed velocity shift 
of 20 m s _1 is comparable with the uncertainty range esti- 




-0.27 -0.24 -0.21 -0.18 -0.15 -0.12 -0.09 -0.06 



Av (km s ') 

Fig. 4. x 2 as a function of the velocity difference Av 
between the Fen A1608 and AA2382,2600 lines for the 
16-, 11-, and 10-component models which are shown, re- 
spectively, in Figs. 1, 2, and 3. The corresponding xlain 
values are equal to 231.631, 135.314, and 166.397. The 
minima of the curves give the most probable values of 
Aa/a (sec eq.[l]): 5.4 ppm (the 16-componcnt model) 
and 4.8 ppm (the 11- and 10-componet models). The la 
confidence level is determined by A% 2 = 1 which gives 
<JAv = 0.080 km s" 1 , 0.065 km s" 1 , and 0.085 km s" 1 or 
(TAq/q = 2.4 ppm, 1.9 ppm, and 2.5 ppm for the 16-, 11-, 
and 10-component models, respectively. 



mated from the ThAr lines. So, in what follows we consider 
A2382 and A2600 as having the same radial velocity, and 
calculate Aw between this velocity and that of the line 
A1608. 

To demonstrate that our result is not sensitive to 
the number of subcomponents, we have analyzed differ- 
ent models. The 16-component model is chosen as one 
which covers the whole Fen profile in the range from 
—90 km s _1 to 55 km s _1 (Fig. 1) and provides xl 55> 1 
with minimum number of components. The other models 
are shown in Figs. 2 and 3 with the corresponding pa- 
rameters listed in Table 2. The 11-component model in 
Fig. 2 ignores the weak absorption in the wings of the 
Fell lines at -90 km s _1 < v < - 55 km s _1 and 40 
km s -1 < v < 55 km s _1 . In the model shown in Fig. 3 
the number of components is reduced to 10 which is found 
to be the lower limit since models with 9 components can- 
not be fitted adequately to the Fen lines in the range —55 
km s^ 1 < v < 40 km s _1 with xl l ess than 3. We note 
that all of the analyzed models show similar radial velocity 
shifts between the A1608 and A2383/2600 lines. The com- 
parison between the theoretical and observed Fe II profiles 
was quantified by the x 2 function (defined by eq.[2] in 
Paper II). To find the most probable value of Av, we fit the 



6 



S. A. Levshakov et al.: Aa/a at z = 1.84 (RN) 



absorption lines with a fixed Av, changing Av in the inter- 
val from —270 m s _1 to —90 m s _1 in steps of 10 m s _1 (see 
Fig. 4). For each Av, the strengths of the sub-components, 
their broadening parameters and relative velocity posi- 
tions were allowed to vary in order to optimize the fit 
and thus minimize x 2 . The \ 2 curve as a function of Av 
in the vicinity of the global minimum must be smooth and 
have a parabola-like shape. The most probable value of Av 
corresponds to the minimum of \ 2 ■ ln our case it is —180 
m s _1 for the 16-component model. The lcr confidence in- 
terval to this value is given by Ax 2 = X 2 — Xmin = 1 ( e -§-; 
Press et al. 1992), which is shown by the horizontal dashed 
line in Fig. 4. It is seen from the figure that the lcr errors 
are equal to 80 m s" 1 , 65 m s _1 , and 85 m s _1 for the 
models in question. Adding quadratically the wavelength 
scale calibration error between the blue and red arms of 
30 m s _1 and using (1), one can easily obtain Aa/a = 
5.4 ± 2.5 ppm for the main 16-component model. 

4. Possible systematics 

We have specifically investigated those effects which might 
introduce a non-zero difference between the blue and the 
red lines and thus simulate a Aa/a variation at the ppm 
level. In addition to the systematics discussed in Sect. 2, 
the following factors should be taken into account. 

Systematic shift between the blue and red arms. 
Different velocity offsets may occur in the blue and red 
frames causing an artificial Doppler shift between the 
A1608 and A2382/2600 lines. To probe this effect we used 
the same procedure as for the QSO exposures and cali- 
brated in wavelength the spectra of the ThAr arcs. For 
the blue and red frames we selected well-exposed ThAr 
emissions which bracket the positions of the Fen A1608 
and A2382/2600 lines. The calibrated ThAr spectra do 
not show the rms errors larger than 20 m s _1 . In Fig. 3, 
an example of the residuals (open circles) between the lab- 
oratory and fitted ThAr wavelengths for the 4th exposure 
is shown. The total number of ThAr lines used to cal- 
culate the rms error in the echelle order is indicated on 
each panel, whereas only the ThAr lines close to the ob- 
served positions of the Fe n lines of interest are depicted. 
A conservative shift between the blue and red arms of 30 
m s _1 can be obtained if the rms error of 20 m s _1 is taken 
for each arm. 

Centering the target onto the slit. Radial velocity shifts 
induced by different changes of isophote onto the slit are 
not chromatic, i.e. the shifts are the same in the blue and 
red UVES arms. The fact that we have re-centered the 
object in the different nights and also before and after 
passage to the meridian assure that the effect was different 
from one exposure to the other, and, hence, averaged out 
in the mean spectrum. Besides, tracking movements act 
in smoothing the isophote onto the slit. We note that in 



our previous results on HE 0515-4414 (Paper II) we also 
used Fe n lines falling in both UVES arms without showing 
evidence of shifts. 

Isotopic shifts. Using the mass shift constants /cms cal- 
culated by Kozlov et al. (2004), we estimated the uncer- 
tainty of Aa/a ~ 0.77 ppm caused by unknown isotope 
abundances (Paper I) which is equivalent to the uncer- 
tainty of about 20 m s _1 . 

Unresolved components. The normalized residuals, 
[Tf al - Ff s )/o-i, in Figs. 1, 2, and 3 do not show series 
which could suggest unresolved Fen components. 

Possible blends. The A1608 and A2382 lines lie in the 
spectral regions free from telluric absorptions. No metal 
absorption lines from other systems which could blend 
with these iron lines were found as well. Weak night sky 
absorptions affect the A2600 profile at v = —30 km s _1 and 
in the range v < —50 km s^ 1 which was not included in the 
fitting of this line. Since both the A2382 and A2600 lines 
are considered as having the same radial velocity and the 
line A2382 is clean, this contamination of A2600 does not 
influence the final value of Aa/a. 

5. Conclusions 

We have re-observed the quasar Q 1101-264 with spec- 
tral resolution FWHM ~ 3.8 km s _1 and re-analyzed the 
Fen profiles associated with the z — 1.84 damped Lya 
system. The data represent one of very few spectra of 
QSOs obtained with this high resolution and S/N > 100. 
The newly obtained value of the relative radial velocity 
shift Av = -180 ±85 m s" 1 between the A1608 and 
AA2382, 2600 lines represents a factor of 1.5 improvement 
with respect to lower resolution archive data (FWHM ~ 
6 km s _1 , Paper I). Admittedly, there might be further 
hidden systematic effects which would challenge the inter- 
pretation as due to variation of a. If real, it would corre- 
spond to Aa/a = 5.4±2.5 ppm. Thus, the highest spectral 
resolution achievable on QSO spectra at 8-10 m telescopes 
is proved to significantly contribute to higher accuracy in 
the Aa/a measurement. 

We note that averaging the 8 systems in the redshift 
interval 1.5 < z < 2.0 of Chand et al. (2004) one obtains 
Aa/a= 2.4 ± 1.2 ppm. On the other hand, Murphy et 
al. (2003) results averaged over the same redshift interval 
(16 systems) give Aa/a = —5.4 ± 4.5 ppm. Whereas the 
Chand et al. and our values are consistent within the lcr 
confidence interval, the value of Murphy et al. is signifi- 
cantly different. 

Up to now, we have obtained high precision measure- 
ments of the Aa/a with the SIDAM method at two red- 
shifts, z = 1.15 and z = 1.84, only. Of course, to probe 
the true redshift evolution of a many more measurements 
are needed. It should be understood, however, that such 
measurements are extremely hard to perform since they 
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Table 2. SIDAM analysis: positions, v (in km s -1 ), of the Fen sub-components, the line broadening parameters, b 
(in km s _1 ), and the column densities N\2 (in the units of 10 12 cm -2 ) for the 16-, 11-, and 10-component models 
shown in Figs. 1, 2, and 3 (the corresponding uncertainties of the fitting parameters are less than 5%). The v values 
are referred to the A1608 A line, whereas the positions of the A2382 and A2600 A lines are shifted at v' = v + Av with 
the optimal values of Av given in Fig. 4. 





V 


b 


N12 


V 


b 


N12 


V 


b 


N12 


No. 


16-component model 


11-component model 


10-component model 


1 


-74.494 


5.36 


0.48 














2 


-68.198 


1.59 


0.11 














3 


-62.067 


2.42 


0.34 
















-39.630 


6.39 


3.92 


-39.643 


6.30 


3.95 


-39.661 


6.30 


3.91 


5 


-37.029 


2.85 


2.36 


-37.018 


2.79 


2.30 


-37.046 


2.82 


2.33 


6 


-24.643 


7.02 


2.90 


-24.650 


7.14 


2.95 


-24.497 


7.30 


3.02 


7 


-17.879 


2.05 


0.99 


-17.837 


2.11 


1.03 


-17.902 


1.86 


0.90 


8 


-8.875 


5.69 


2.02 


-9.677 


4.67 


1.58 


-8.667 


6.13 


2.11 


9 


-1.279 


3.60 


9.56 


-1.171 


3.88 


10.60 


-1.072 


3.76 


10.15 


10 


4.085 


3.31 


2.67 


4.607 


3.03 


2.00 


4.639 


2.92 


2.01 


11 


15.493 


2.77 


0.29 


15.504 


3.02 


0.31 


16.905 


4.47 


0.47 


12 


20.110 


1.18 


0.28 


20.352 


1.50 


0.34 








13 


25.089 


3.03 


2.31 


25.145 


2.89 


2.23 


24.873 


3.45 


2.43 


14 


35.192 


2.95 


2.25 


35.232 


3.10 


2.31 


35.288 


2.95 


2.27 


15 


42.746 


2.36 


0.33 














16 


47.827 


2.31 


0.48 















require spectral data of highest quality which can be ob- 
tained for the brightest QSOs only. In spite of the fact 
that many giant telescopes were put into operation in the 
last decades, the unique simultaneous wavelength cover- 
age 3000 A < A < 10000 A of UVES makes it the only in- 
strument suitable for this work. Thus, the progress in our 
understanding the evolution of Aa/a and, hence, of dark 
energy is expected to be 'slow- rolling' until more powerful 
Extremely Large Telescopes will come into operation. 
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Fig. 5. An example of the residuals (open circles) be- 
tween the laboratory and fitted ThAr wavelengths for the 
4th exposure. The total number of the well-defined ThAr 
lines used to calculate the rms error in the echelle order 
is indicated on each panel, whereas only the ThAr lines 
close to the observed positions (marked by vertical bars) 
of the Fen lines of interest are depicted. 



